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SPECTRORADIOMETRIC  ANALYSIS  OF  RADIO 

SIGNALS. 

By  Chester  Snow. 


ABSTRACT. 

When  a  radiation,  regularly  repeated,  is  received  by  a  simple  circuit  possessing 
inductance,  capacity,  and  resistance  there  is  induced  in  it  an  electric  current  whose 
effective  value  may  be  measured  by  a  hot  wire  ammeter,  provided  the  signal  is  re- 
peated for  a  sufficient  time  to  take  readings  of  the  instrument,  and  provided  that  the 
time  between  the  signals  is  very  small  compared  to  the  period  of  the  ammeter. 
If,  while  the  radiation  continues,  the  circuit  be  tuned  to  various  frequencies  in  a  cer- 
tain range,  and  the  current  observed  in  each  case,  then  these  current  readings  may  be 
plotted  against  resonance  frequency  of  the  circuit.  If  the  latter  has  a  relatively 
high  resolving  power  compared  to  the  spectral  structure  of  the  wave,  such  a  curve 
will  be  proportional  to  the  spectral  distribution  of  energy  in  the  radiation.  As  a 
matter  of  fact  this  is  not  usually  the  case,  and  to  derive  the  energy  distribution  of 
the  wave  requires  that  a  certain  integral  equation  be  solved  in  which  the  above- 
mentioned  curve  is  the  known  function.  Several  methods  of  obtaining  a  practicable 
solution  of  this  equation  are  developed  here  and  applied  to  experimental  observations. 
The  function  of  this  paper  is  to  develop  a  method  of  spectrum  analysis  which  may  be 
applied  to  the  study  of  the  interference-producing  quality  of  a  given  sending  station. 
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I.  INTRODUCTION— DERIVATION   OF   FUNDAMENTAL 

EQUATION. 

The  question  as  to  what  constitutes  interference  quality  in  a 
radio  wave  is  one  of  such  importance  as  to  deserve  both  experi- 
mental and  theoretical  study.  Undoubtedly  a  thoroughgoing 
spectrum  analysis  of  the  radiation  when  properly  interpreted 
ought  to  furnish  a  logical  basis  for  rating  the  interference- producing 
qualities  of  a  sending  station.     It  must  be  confessed  that,  in  spite 

231 


232  Scientific  Papers  of  the  Bureau  of  Standards  [Voi.  iq 

of  the  urgency  of  the  need  of  it,  no  satisfactory  meaning  for  the 
term  "coefficient  of  interference  of  a  sending  station"  has  yet 
been  arrived  at,  either  from  a  scientific  or  a  legal  aspect.  The 
accepted  notion  of  "decrement"  seems  to  have  been  considerably 
overworked,  and  it  is  hoped  that  the  somewhat  larger  conception 
of  tliis  article  may  prove  useful. 

Some  time  ago  R.  T.  Cox,  of  this  bureau,  began  an  experimental 
analysis  of  certain  simple  types  of  radiation  with  the  object  of  devel- 
oping a  method  by  which  the  interference-producing  quality  of  the 
wave  might  be  characterized.  In  interpreting  the  results  obtained, 
it  immediately  became  evident  that,  even  with  a  circuit  of  the 
finest  attainable  resolving  power,  the  record  partakes  of  the 
nature  of  the  receiving  circuit  as  well  as  of  that  of  the  incoming 
wave.  To  disentangle  the  latter  from  the  former  requires  the 
solution  of  a  certain  integral  equation,  and  this  constitutes  the 
aim  of  this  article. 

If  the  radiation  from  a  given  sending  radio  station  is  received 
simultaneously  by  a  set  of  simple  circuits  whose  constants  are 
properly  chosen,  the  arrangement  constitutes  a  spectroscope. 
From  measurements  of  the  effective  current  induced  in  each  cir- 
cuit, the  spectral  distribution  of  energy  in  the  wave  may  be  found. 
Or  if  the  message  can  be  repeated,  a  single  circuit  may  be  used 
and  its  constants  varied,  the  effective  square  of  the  current  being 
measured  for  each  value  of  the  constant  which  characterizes  the 
circuit. 

When  the  circuit  possesses  very  high  selectivity  so  that  it  may 
be  considered  to  respond  to  a  definite  frequency  only,  then  if  it 
be  tuned  successively  to  various  frequencies  and  the  mean  square 
of  the  current  be  measured  for  each,  this  will  be  proportional  to 
the  spectral  distribution  of  intensity  in  the  radiation.  Such  an 
ideal  case  is  not  always  attainable  in  practice.  The  observed 
current  as  a  function  of  frequency  is  often  profoundly  affected  by 
the  character  of  the  receiving  circuit.  The  resonance  curve  of 
every  circuit  possess  a  finite  width,  due  to  the  fact  that  there  is  a 
dissipation  of  energy  in  the  circuit.  A  spectral  structure  in  the 
radiation,  which  is  fine  compared  to  the  resonance  curve  of  the 
circuit,  may  be  entirely  masked.  The  problem  which  then  pre- 
sents itself  is  to  find  how  the  circuit  impresses  its  character  upon 
the  record  in  order  to  make  corrections  for  it,  and  draw  inferences 
as  to  the  character  of  the  radiation  which  may  be,  at  first  glance, 
entirely  hidden  in  the  observations.     It  is  analogous  to  attempt- 
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ing  to  resolve  a  light  source  b}'-  means  of  a  spectroscope  of  insuffi- 
cient resolving  power,  where  one  makes  up  for  lack  of  resolving 
power  by  a  quantitative  study  of  the  energy  distribution  in  the 
record,  a  procedure  which  is  feasible  in  the  range  of  radio 
frequencies. 

The  current  in  a  circuit  possessing  resistance,  inductance,  and 
capacity,  R,  L,  and  C,  respectively,  oscillating  freely,  is  of  the 
form  Ae'^'^^^  sin  27rx  (t  +  y) 
Where 


R         ,         1        I       R'  .  . 

a  =  — irandaj  = — -1/77=. ri  (i) 

4xL  27r  V  LC     4U  ^  ' 

The  quantity  5  = is  usually  called  the  logarithmic  decre- 
ment of  the  circuit.  If  the  capacity  only  is  varied,  the  damping 
coefficient,  27ra,  remains  constant,  and  the  characteristic  fre- 
quency, X,  may  be  given  any  value,  real  or  imaginary.  It  will 
become  evident  that  if  a  radiation  is  examined  for  all  real  values  of 
X,  this  is  sufficient  to  determine  its  energy  distribution. 

The  wave  produces  in  the  circuit  an  electromotive  force  (g(i)  = 
En  sin  2x  %n  (^+o;n)  where  £^n  is  referred  to  in  the  following 
as  the  Fourier  component  of  the  radiation.  It  is  proportional 
to  the  electromagnetic  field  in  the  case  where  the  receiving  circuit 
consists  of  antenna  with  capacitative  coupling  to  earth,  but  in  the 
case  of  coil  aerials  E^  will  be  proportional  to  the  time  derivative 
of  the  field  and  hence  contains  the  frequency  as  a  factor.  This 
will  maintain  in  the  circuit  the  current  3  (0  =/n  sin  2  -k  x^  (^  +  |8n) 
whose  mean  square  is 

/  2     y  17  2 

^    ~    2   ~2 


\  27rXnC/ 

F  2 
8^L'[4a'+(.„-fJ] 


where  X  =  — -./^r— =  the  resonance  frequency  of   the  circuit,  or 

the  frequency  which  would  elicit  the  greatest  response  from  the 
cii'cuit.  It  is  connected  with  the  characteristic  frequency,  x,  by  the 
relation 

a;2  =  X^-a2  (2) 
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There  is  very  little  difference  between  these  two  frequencies 
for  a  good  receiving  circuit,  but  the  exact  relations  which  follow 
are  more  simple  when  expressed  in  terms  of  x.     The  frequency, 
X,  is  the  one  given  by  wave-meter  measurements. 
If  for  brevity  we  let 

/  (x)  =SRL  ixW')  (3) 

then 

^ax  E^a 


/w  = 


""    4«^  + 


/        a^  +x^V 


for  a  simple  harmonic  wave  of  frequency  Xn  and  amplitude  (of 
emf)  En. 

A  similar  formula  holds,  however,  if  the  electromotive  force  is 
expressed  as  a  finite  or  infinite  number  of  such  components,  if  3^ 
denotes  the  time  average  of  current  squared  over  a  large  number 
of  periods.  Making  the  usual  assumption  that  it  holds  when 
there  are  present  all  frequencies  from  x'  =0  to  ic'  =  «  gives  the 
integral  equation 

oo 

,        4ax  r       E^  (x')  dx' 

^  ^""^  ~  "^  J  TJ77J~^^+^'  (4) 


(^'-"-^7 


'6  4^'  + 

to  determine  E^  (x)  in  terms  of  the  function  /  (x)  which  is  given  by 
observations. 

The  proof  in  the  case  of  a  Fourier's  series 

3  W  =  2  /„  sin  ^  (t  +/3„) 
«=i  •' 

is  effected  by  squaring  the  series  and  noting  that  the  integral 

NT 

j^  I  sin  -^  {t  +j8„,)  sin  -y-  (t  +^n)  dt 

0 

taken  over  a  large  number  of  periods,  NT,  is  practically  zero  if  m 
is  not  equal  to  n  and  is  one-half  if  m  =n.     This  gives 

I\     I  ^  E\ 


NT 
0 


3'=Krr        3=*  Wrf/=S-f  =  -S 


R'+llTTXnL ^.Y 

\  27rXnC/ 
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and  holds  no  matter  how  large  the  finite  value  of  T  which  measures 
the  range  of  convergence  of  the  series.  It  is  assumed  to  hold  even 
when  the  range  of  convergence  is  infinite  and  the  series  passes  into 
a  Fourier  integral. 

The  integral  in  equation  (4)  involves  x^  as  a  parameter  under 

the  integral  sign,  so  that  -  f(x)  is  really  a  function  of  x^  (and  of  a). 

Therefore,  E'^  will  be  determined  as  a  function  of  x'^,  and  E^ 
{—x')  =E^{x').  Since  E'^{x')  is  an  even  function  of  x'  the  equa- 
tion (4)  may  be  written 


/w  = 


E\x')dx'  \       x'E^{x')dx'^ 


—  I  /  2+  2V  =  —  I  


00  00 

^a_  rx'EKx')dx'       a   Cx'E'{x')dx' 

2irJ  a^+{x'-xy     2TrJ  a'+i^x'  +xy 


If  the  variable  of  integration  be  changed  in  the  second  integral, 
letting  x'  =  —  IX,  it  becomes 


+  - 
27r 


a_  rnE'(-fi)dfi_  a   C x'E^{x')dx' 
2TrJ  a^  +(fi  —  x^)  ~ 2irj  a^  +{x'  —  x)^ 

—00  —00 

Hence,  if  we  let  u(x)  =x  E^(x)  the  fundamental  equation  becomes 

CO 

.        a  r    u{x')dx' 

'^^^~Tr}  a^+{x'-xy  ^5; 

—00 

where  it  is  known  that  fix)  and  u{x)  are  odd  functions  of  x.  The 
equation  (5)  formulates  the  problem  here  undertaken,  which  is  to 
find  the  unknown  function  u{x)  in  terms  of  the  observed  or  given 
function  f{x).  The  remainder  of  this  paper  consists  in  the  de- 
velopment of  several  forms  of  solutions  of  this  equation,  together 
with  illustrative  examples  and  applications  to  experimental 
results. 
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11.  FIRST    FORM    OF    SOLUTION— FOURIER'S    INTEGRAL 

METHOD. 

It  may  be  verified  by  direct  integration  that 


>""' 


cos  oc  {%'  —  x)  da  ii  a  >  o 


a^  +  (x'-xy 

0 

Hence,  the  equation  (5)  may  be  written 

-oo  00 

f{x)=-  \u  (x')  dx'  le-'"'  cos  a  (x'-x)  da  (6) 

—00  0 

Comparing  this  with  the  Fourier  integral  identities  for  /  (x)  and 
u  (x) 


00  00 

/  (^)  =i  Cj  (x')  dx'  fcos  a  (x'-x)  da 


(7) 


u  (x)  =-  \  u  (x')  dx'  j  cos  a  {x'  —  x)  da  (8) 

leads  to  the  inference  that 


00  00 

u  {x)=-\f  {x')  dx'  j  e^"  cos  a  (x'  - x)  da 


(9) 


or 


oo  CO 

u  (x)  =-  \  f  (m)  dfj.  I  e'^"  sin  afj.  sin  axda  (10) 

0  0 

The  latter  form  is  a  consequence  of  the  fact  that  f(x)  is  an  odd 
function  of  x.  That  the  formula  (10)  does  give  the  solution  of  (5) 
may  be  verified  by  making  use  of  Laplace's  integrals  ^ 

00  00 

J  cos  a  tj,       Ct  sin  a.  t,.  -,„  -r         ^  /     n 

—  00  —CO 

Replacing  x  by  x'  in  (10)  and  substituting  in  the  right  hand  side 
of  (5)  gives  1 


CO  03  CO 

2a  fx  /  N  J      r  i»    •  J      fsin  a  x'  dx' 


2 

? 
"0 


'  Meyer,  Bestimmte  integrale,  p.  19 
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To  perform  the  integration  first  with  respect  to  x'  let  x'  -x  =t. 

Then 

00  00 

J  a^+{x'-xy         J      a"  +f 

—CD  —00 

00 


,  .  COS  a  t,.     T  e^"    . 

=sm  a  x  I  -^^ — 5 — -^dt  = sm  a  x 

+f  a 


The  triple  integral  last  written  then  reduces  to 

00  oo 

-  I  /  (^i)  £^  /i  I  sin  a  n  sin  a  x  d  a 
0  0 

which  is  identically  equal  to  /(x)  and  therefore  (10)  is  the  solution 
of  (5).  The  understanding  in  (10)  is  that  the  integration  is  to  be 
performed  first  with  respect  to  x',  otherwise  it  is  meaningless. 
This  form  of  solution  has  the  advantage  that  it  is  applicable  when 
f{x)  has  any  finite  number  of  finite  discontinuities.  When  the 
observed  function  f{x)  is  given  by  a  curve,  the  formal  solution  (9) 
is  unsatisfactory.  It  is  desirable  to  work  out  a  method  whereby 
the  function  u(x)  may  be  derived  from  the  curve  for  f{x)  by  the 
use  of  graphical  or  mechanical  methods.  This  is  done  in  the 
following. 

III.  SECOND     FORM     OF     SOLUTION— ADAPTATION     FOR 
EVALUATION     WITH  A     HARMONIC     ANALYZER. 

In  practice  we  are  usually  concerned  with  a  relatively  small 
range  of  high  frequencies,  to  which  the  observations  are  confined. 
As  far  as  these  observations  are  concerned,  we  may  assume  that 
u(x')  =0  outside  this  range  x^Kx'  <  x^.  It  is  to  be  noticed  that  in 
the  original  equation  (5)  the  nucleus  involves  {x  —  x'),  the  differ- 
ence of  frequencies,  and  therefore  we  may  take  the  origin  of  x 
anywhere,  and  the  form  of  (5)  and  of  the  solution  (10)  will  be 
unaltered.  Let  us  take  it  at  the  middle  of  the  observation  range, 
which  is  supposed  to  extend  a  distance  ±  T  from  this  origin.  The 
curve  f{x)  having  been  plotted  according  to  its  definition  ]{x)  = 
S  R  L  xQ"  where  x  is  actual  frequency,  we  shall  denote  it  by 
F(x),  where  x  now  represents  frequency  difference  from  the  new 
origin.  Similarly,  we  shall  use  v{s)  instead  of  u(x')  for  the  func- 
tion which  determines  the  energy  of  the  wave  as  a  function  of 
this  frequency  difference. 

58954°— 23 2 
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Since  v{s)  is  zero,  if  ^<  —  T  or  ^>T,  the  fundamental  equation 
to  determine  energy  -v^s)  when  the  observed  curve  Fix)  has  been 
plotted  is 

T 

„,  ,      a  C     v{s)  ds  I     V 

ttJ  a^  -\-{s  —  %y 


-T 


where  v{x)  and  F{x)  are  no  longer  known  to  be  either  even  or 
odd  functions  of  x.  The  assumption  that  v(s)  is  zero  outside  the 
range  x^Kx'  <  X2  is  contradictory  to  the  fact  already  noted  that 
E^{x')  is  an  even  function  of  x'.  In  keeping  with  this  fact  there 
should  be  added  to  the  right  side  of  equation  (12)  the  corre- 
sponding "  image  "  of  the  term  there  present  corresponding  to 


IT  J  a"" 

—  X2 


u{x')  dx' 
+  {x'-xy 

which  may  be  put  in  the  form 

){s)  ds 
+  {s+xY 


aC    V 

~T 


which  is  utterly  negligible  when  x  is  large  and  positive  compared 
to  the  term  retained  in  (12).  The  nucleus  of  the  integral  in  (12) 
may  be  written  K{x  x')  where 


IT  a''  +{x  —  yy 

This  has  a  sharp  maximum  —  at  :*;  =y,  reduces  to  half  this  value 

^  ira  -^ 

when  x  —  y=a  (the  "half  width  of  the  resonance  peak")  and 
rapidly  decreases  when  (x—y)  becomes  large  compared  to  a. 
This  shows  that  the  fraction 


ira^+ix+yy 


must  be  negligible.  If  the  wave  is  undamped,  its  spectrum  is  a 
single  line  at  the  frequency  x^  say.  The  corresponding  observed 
ftmction  is 

a         A 

Ta'+ix—x^)' 
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A 
which  has  the  maximum  value  —  at   the  resonance  frequency 

ira 

x-=Xj^  and  reduces  to  half  this  value  when  x=Xi±a.  Hence  a 
is  called  the  half  width  of  the  resonance  curve.  If  the  spectrum 
consists  of  two  equal  lines  separated  by  the  distance  Ax,  in 
frequency,  the  observed  curve  would  be  the  sum  of  two  such  terms 

„,  .     A  a  A 


T  a^  +(x  —  x^)^     Tra^+ix  —  Xi)^ 

and  if  |a;i  — iC2|>  a  this  curve  would  have  two  visible  peaks,  which 
would  coalesce  (practically)  into  one  peak  if  Ax  or  \x2—Xi\^a. 
Therefore,   a  measiu^es  the  limit  of  resolution  of  the  receiving 

circuit  as  an  analyzer  and  —  =  -  =  —  is  the  resolving  power  of 

■^  X  X       2T 

the  circuit  as  a  spectroscope.  In  the  notation  immediately 
following  the  limits  will  not  be  written,  but  it  is  understood  that 
the  integrals  are  all  definite  integrals  over  the  range  from  s  =  —T 
tos  =  +T. 

The  general  theory  of  integral  equations  shows  that  the  solution 
of  the  integral  equation 

F{x)=fK{xs)v{s)ds  (13) 

may  be  found  in  terms  of  the  normal  functions  4>n  ix),  and  the 
associated  functions  ;^n  (x)  and  the  characterisitc  constants  Xn 
where 

y*^n  W  0m  {x)  dx  =0  if  ny^m,  =1  if  n  =m  (14) 

and  <^n  (x)  and  r^n  (x)  are  solutions  of  the  pair  of  simultaneous 
homogeneous  integral  equations. 

fn  (x)  -\nfK  (xy)  0n  (y)  dy  =0 

(15) 

0n  (x)  -\nfK  (y  x)  xpn  (y)  dy  =0 

and  where  the  constants  Xn  are  the  roots  of  the  characteristic 
determinant  D  (X)  =0  which  is  found  when  K  (x  y)  is  known  b}^ 
Fredholm's  formula. 

If  the  nucleus  K  {xy)  is  a  symmetric  function  of  x  and  y,  so  that 
K  {xy)  =K  iy  x) ,  then  i/'n  {x)  becomes  identical  with  0n  {x) . 

If  the  roots  of  D  (X)  =0  are  infinite  in  number,  there  are  an 
infinite  number  of  normal  f  mictions  ^n  (x)  and  the  associated  ones 
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rpn  (x) .  In  this  case  any  function  v  (x)  may  be  developed  in  a  series 
of  these  normal  functions,  analogous  to  Fourier's  series. 

V   (x)   ='EVn  </>n   (x)  (16) 

provided  "Ev^a  is  convergent.  The  constant  coefficient  Va  is  found 
by  multiplying  this  equation  by  4>n  (x)  dx  and  integrating  over  the 
complete  range.     This  gives 

T^n  =fv   (X)    0n    (X)    dx.  (17) 

Multiplying  the  first  of  the  equations  (15)  by  ^^  (x)  dx  and  integrat 
ing  gives 

y*  }pn  (x)  ^m  (x)  dx=\af  (f>a  (y)  dy  f  K  {x  y)  ipm  (x)  dx 
=f  <^n  iy)  <^m  (y)  dy  by  (15) 

which  shows  that  the  i/'n  (x)  themselves  constitute  a  set  of  normal 
functions.     Consequently 

F  (x)  =S  F„  lAn  (x)  where  F^=^fF  {x)  ^^n  {x)  dx  (18) 

n 

The  work  of  Hilbert  and  Schmidt  has  shown  that  the  nucleus 
may  be  expressed  in  terms  of  these  functions  by  the  formula 

K{xy^-E'tM^  (19) 

To  solve  the  original  integral  equation,  assume  that  the  un- 
known function,  which  is  to  be  found,  may  be  expressed  by  a 
convergent  series  of  the  form 

v{x)  =2  i;k0k(aj)  (20) 

k 

The  solution  consists  in  finding  the  coefficients  v^.  To  do  this 
substitute  this  series  for  v{x)  in  the  integral  equation,  substituting 
at  the  same  time  the  series  for  F{x)  in  terms  of  the  associated 
functions  t/'n(x).     Equation  (13)  becomes 

F{x)  =SF„^n(^)=S^S'z;k  r«^„W0kWrf^  =  2^VnW 

n  n      ^n      k        J  n'^n 

This  equation  will  be  satisfied  if  r-^  =  Fn  or 

Vu  =XnF„  ^KfF(s)Ms)ds  (by  18) 
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Hence,  the  solution  is 

v{x)  =2  KFncf>n{x)  =  S  X„0„(x)  CF{s)Ms)  ds  (21) 

n  J 

provided  this  series  is  convergent.  The  coefficients  F^  must 
first  be  found  and  then  the  sum  of  the  series  taken.  An  exceed- 
ingly great  gain  in  convenience  would  be  made  if  it  were  possible 
to  interchange  the  order  of  limits  in  (21)  so  that  we  could  sum  the 
integrand  first  and  then  integrate  with  respect  to  x  giving 

V  {x)  =  fF  is)  ds  S  Xn  0„  (x)  i/'n  (s)=fF  (s)  N  (x  s)  ds 

(22) 
where 

A^  (xy)=ZK(t>,  (x)  xl^n  (y) 

for  this  function  N  could  be  computed  once  for  all  for  a  given 
receiving  circuit  and  a  given  range  of  integration,  and  the  work 
of  applying  it  to  various  observed  functions  F  (s)  by  (22)  would 
involve  considerably  less  routine  work  than  in  applying  (21). 
However,  it  is  to  be  noticed  that  this  series  for  A^  {x  y)  can  sel- 
dom converge,  for  Xn  increases  with  order,  n,  and  if 

2   0n    (X)    ^n    iy) 
n  Xn 

converges  it  is  not  to  be  expected  that  S  Xn  0n  (x)  \pn  (y)  will  con- 
verge in  general. 

In  order  that  the  preceding  developments  shall  be  valid,  it  is 
necessary  that  the  system  of  normal  functions  be  "closed,"  which 
will  be  the  case  if  no  function,  h  (x) ,  exists  such  that  X  K  {x  s)  h  (s) 
ds  =  o.     The  normal  fimctions  are  then  infinite  in  number. 

The  above  form  of  solution  is  of  no  practical  use  for  computing 
unless  one  can  find  the  set  of  constants  W  ....  and  the  correspond- 
ing normal  functions. 

However,  in  case  the  nucleus  is  a  function  oi  x  —  y  or  oi  x  +  y,  or 
is  the  sum  of  such  functions,  the  normal  functions  are  trigonomet- 
ric functions.  The  harmonic  analysis  of  the  nucleus  (either 
mechanically  or  analytically)  then  amounts  to  a  determination  of 
the  characteristic  constants.  In  the  present  case  K  (x  y)  is  an 
even  function  of  {x  —  y)  and  may  be  resolved  as  follows : 

K{xy)  =  '-  i^o  +S  i^n  cos  VlILS^ 

2  re  1 
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or  if 


T 


.  n  IT  u   J 

ttT 


d'+iJi' 

-T 


then 


K  («>/)=- /Co +2  Kn  cos  ^ 


n  T  X         n  TT  y 
cos  —~- 


,„„      .     n  TT  X    .     n  IT  y 

+  S  An  Sin     ^     sm      „-^ 

«  ^  ■* 

Also  F  (x)  may  be  developed  in  a  Fourier's  series  .... 
T<x<T 

F  {x)=-  Fo  +2  Fn  cos  -^jT-  +2  F„'  sm  -^^r- 
2  „  i  i 

where 


and 


^°  ""  T  I  ^  ^'"^  ^^^  ^^^^  ^  '^ 


(23) 


K24) 


F„'  =  :i  1  F  (m)  sin  ^^^^  d  ju 

-T 

If  a  value  of  v{x)  is  sought  in  the  form 

/  \      I        vn  ^Trrt     _,    ,     .    kirx  ,     . 

DXx)  =-^0  V^Vk  cos  -7fr  +  27;&'  sm  -7fr-  (25) 

2        fe  i        k  i 

if 

-r<:r<r 

the  coefficients  Vo,  Vn,  and  i;n'  may  be  found  by  substituting  all  these 
series  in  the  equation 

T 

■^W  =fK{x  y)v{y)  dy 

-T 


which  gives  the  equation 


-Fo  +2F„  cos  -™-  +  2F'  sm  -„- 

: 7  +7^  2iC„'z;n  cos  ^^  +  T  2/Cn^'  n  sm  ^=r- 

2  „  1  „  I 
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This  will  be  satisfied  if  we  take 

"Vo ■'^  o  __    "  n  /     __  ^    n 

'^"^Kof         ^"""KJ-         ^°~KnT 
and  the  solution  is 

I  fi  Fo  ,  ^Fn         nirx       F„'    .    nirxl 

if  this  series  is  convergent. 

The  formula  (26)  gives  a  practical  method  of  solving  the 
problem.  The  coefficients  Kn  may  be  evaluated  once  for  all  for 
a  given  circuit  (a  given)  and  a  given  range  2  T.  However,  this 
must  be  done  by  approximate  methods,  or  mechanically,  since 
the  integral  which  defines  Kn  can  not  be  integrated  in  finite 
terms  if  T  is  finite.  While  it  would  be  very  convenient  to  reverse 
the  order  of  limits  in  (26)  and  write 

■v{x)=fF{s)  N  (x-s)  ds 

-T 

where 

nirx 
_cos-=- 


it  is  easy  to  see  that  such  a  procedure  is,  in  general,  not  allowable, 
since  this  series  will  not  converge  in  general.  It  is  easy  to  show 
by  means  of  Laplace's  integrals  that  the  limit  of 

KnT  =  ire    ^  when  F  =  » 

and  if  we  let  dx  =  -jfic  =  -^  the  solution  (26)  approaches  as  a  limit 
when  F  =  00 

00  cx>  . 

v{x)  =  -  I  F{s)ds  I  e^"  cos  a{s  —  x)da 
00  0 

which  has  a  meaning  only  if  the  integration  is  performed  first 
with  respect  to  s.  This  is  identical  with  the  solution  previously 
obtained  (9) 
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IV.  THIRD  FORM  OF  SOLUTION— OPERATION  WITH 

cos  a  D^. 

By  multiplying  the  identity 


00 


a 


g-aa  (,Qg  ci  (x'  —x)  da  =- 


'ira'+{x'-xy 
by  f{x')  dx'  and  integrating  one  obtains  the  identity 

CO  c»  00 

^ J/  {x')  dx'je-^"  cos  a(x'-x)da  -^ J^f^^^T^^  =0  (27) 

—  CO  0  —00 

Adding  this  identically  zero  quantity  to  the  right  side  of  (9)  gives 

CX)  03 

u{x)=-  I  f(x')  dx'  I  (e^"  +e-^")  cos  a  (x'-x)  da 

\  (28) 

_a  r    fix')  dx' 
ttJ  a'+(x'-xy 


The  last  integral  in  this  expression  may  be  evaluated  graphi- 
cally or  mechanically.  It  is  desirable  to  find  a  method  of  evalu- 
ating the  first  integral.     For  this  purpose  it  may  be  noticed  that 

(-  lya'^D:^-"  (a  a)'° 

^^ cos  a  X  =— - — —  cos  a  x 

I 2n  |2W 

and  hence  if  we  symbolize  the  infinite  series  of  differential  opera- 
tions 

I r-^  +  — -^  ....  or  S  ^^ -, ^  by  cos  a  D^    (29) 

[2_  U  n=o  |2W 

then 

2  cos  a  Dj  =  (e'^^'  +  e'^^') , 
and 

2  cos  a  D^  cos  ct  (:c-ic')  ={e^"  +e'^")  cos  a:  (aj-x') 

With  this  in  mind,  the  form  of  the  first  integral  in  (28)  suggests 
that  it  is  the  result  of  applying  the  above  operator  to  the  Fourier 
double  integral  for  f{x),  equation  (7),  applying  the  infinite  series  of 
differentiations  with  respect  to  x  inside  the  sign  of  integration,  a 
process  not  always  allowable    even  with    differentiations  of    a 
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finite  order.  Nevertheless  this  suggests  the  second  important 
form  of  solution. 

CO 

U{X)  =  2  COS    a  D,f  (x)  -  -J  Jl^!^r_^y  (3©) 

—  CO 

This  form,  if  coiTect,  may  often  be  more  convenient  to  evaluate 
than  the  first  form  of  solution  (9),  since  the  operations  on  the 
right  of  (30)  may  all  be  performed  mechanically.  By  applying 
Leibnitz's  theorem  for  the  derivative  of  a  product,  one  finds  that 

T^  X  ^\x)  =x  r^Q'  (x)+-. 3^  (x) 

\2n      ^  "•  ■'         \2n       ^  ■^      \2n-1        ^  ' 

so  that 

-  cos  aD^x'^''  {x)  =2  cos  a  D^  3^  {x)  -  —  sin  a  D^  3^  {x) 

Also  it  may  be  noted  that  for  large  positive  values  of  x'  and  x 
the  term 


a"  +{x'  -\-xy 

is  utterly  negligible  compared  to 

a 
a'+(x'-xy 

and,  therefore,  for  any  actual  range  of  frequencies  encountered 
in  practice  the  energy  distribution  will  be  given  by 

£■2  (x)  20. 

-g^  =  2  cos  a  D^  32  (^)  _  _  sin  a  D^  32  (^) 

CO 

_a^  Cx'  3^  {%')  dx'  (31) 

irxj   a'+ix'-xY 


or 


2ar 


^hD.3M^)-^Vw+. 


(32) 


CO 

a^  Cx'  3"  {x')  dx' 
TTXJ   a'+(x'-xy 


58954°— 23 3 
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provided,  of  course,  that  (30)  is  correct.     In  some  cases  one  or 
two  terms  of  the  series  are  sufficient. 

Having  arrived  by  intuitive  reasoning  at  a  form  of  operation 
which  will  simplify  the  solution  of  the  original  equation,  we  may 
now  inquire  more  carefully  as  to  the  result  of  applying  the  operator 
2  cos  h  D^  to  the  original  equation  (5)  where  o<6<a.     Since 

,  =  real  part  of 


X  a^  +{x' -xY 
the  equation  (5)  may  be  written  in  the  form 


fix)  =  real  part  of  — r  I  -7^ — - — :- 
'    '  ^  TTiJ  X  —x  —  ta 


Assume  that  at  the  point  x,  2  cos  b  D^^  f(x)  has  a  meaning;  that 
is,  the  result  of  applying  the  real  operator  2  cos  b  D^  to  f{x)  gives 
a  single  valued  finite  result,  though  not  necessarily  a  continuous 
function  of  x  and  b. 

Regarding  f(x)  as  defined  by  the  integral  in  the  above  equation, 
2  cos  b  D^  f(x)  is  the  value  of  2  cos  b  D^  applied  to  the  integral, 
and  this  operation  may  be  performed  inside  the  integral  sign 
when  it  leads  to  finite  integrand  and  a  uniformly  convergent 
integral.     Since 

cos  6  I>x  or  S  ^^ '- ^— 

n=o  I  2« 

is  a  real  operator,  the  result  of  applying  it  to  the  real  part  of  a  com- 
plex function  must  be  equal  to  the  real  part  of  the  result  of  the 
operation  upon  the  complex  function.     It  is  easily  seen  that 


|2W 


/        I        \_      (-i)°6^° 
\x'  —x  —  ia)     (x'—x  —  iay' 


Hence, 


2         .  ^           I  2    -  (iby- 

— .  cos  0  Djr  —, r-  =  — .    2 


Tti  ^  x'  —X  —  ia     iri „=o(x'  —x- ia) ^"+* 


provided  this  series  is  convergent.  The  series  will  converge  if  b 
is  less  than  the  modulus  of  {x'  —x  —  ia) ;  that  is,  -^ja^  +{x'  —x)^  and 
the  least  value  of  this  is  a  (when  x=x').     Hence,  if  6<a,  the 
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series  will  converge  for  every  positive  value  of  x  and  x'.     It  is  a 
geometrical  series  and  is  readily  seen  to  be  the  expansion  of 


2 
wi 


X'  —x  —  ia 


.(    '''  .Y 

\x'  —x  —  ta/  _ 


or 


I  r  I  I  1 

Tt\_x'  —x—i  (a  —  b)     x'  —x~i  {a  +b)_] 


the  real  part  of  which  is 
a  —  b  I 


1  + 


a  +  b 


X   {a-by+{x'~xy'    T  {a+by  +  {x'-xy 

consequently 

00 

hF)  fM      ^^~^)  f  u  jx')  dx' 


2  cos 


•*) 


{3Z) 


+  ■ 


a  +  b 


±b  r        u{x')dx' 
T    J  {a  +  by  +{x'-x. 


Now  it  is  possible  to  prove  by  direct  integration  (and  it  will  be 
proven  later  as  an  easy  special  case)  that 


CO 

ab  C  _ 


ds 


a  +  b 


+  {x'-syw  +  {s-xy]     TT   {a+by  +  {x'-xy 


(34) 


Hence,  if  we  change  x  to  ^  in  (5)  and  multiply  the  equation  by 

b  ds 

Tb'  +  (s-xy 
and  integrate,  we  get 

CO  00  00 

b  r    f(s)ds     _^ab  r        ds  r    u{x')dx' 

ttJ  b\  +  s-xy     w'J  b'  +  is-xyj  a'  +  (x'-s 

—00  —CO  —00 

00  00 

-^j  u(x')dx'j  i^h-^+^s-xyUa^  +  ix'-sy] 

—00  —00 

CO 

^a  +  b  r        u  (xQ  dx' 

T  J  {a+by-\-{x'~xy  <-^y34) 


sy 
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and  this  is  the  last  integral  on  the  right  side  of  (33).     Its  value  is 
therefore 


00 


{s)  ds 

+  {s-xy 

or  changing  the  name  of  the  variable  of  integration 

00 
b_  C    f{x')dx' 

ttJ  h'  +  ix'-xY 

—00 

Therefore,  the  result  of  operating  upon  (5)  by  2   cos  hDx  where 
o<&<  a  gives  by  means  of  (33) 

00  00 

{a-h)  C         u{x')dx'  ;  n  ^z  \      ^  C     f{^')dx'        ,     . 

^^ / ,,0       ,    , r-,=  2    COS    hUJix) ■„  '^ ,   ' r^       (35) 

TT     J  {a-hY  +  {x' -xY  '^  ^      7rjb^  +  {x'-xY 

—00  —CO 

By  letting  x'—x={a  —  b)   t  in  the  left  side  of  this   equation,  it 
becomes 


I  ru[x  +  (a 


zMm 


and  by  taking  h  sufficiently  close  in  value  to  a,  this  may  be  made 

u          1      1                   1          i    u{x +0) +u{x  -  o)  ., 

to  approach  as  closely  as  we  please  to  — ^^ ^^ or  u{x) 

if  u  is  continous  at  x. 

The  integral  on  the  right  side  of  (35)  is  easily  proven  to  be  a 
continous  function  of  h  and  x  even  if  f{x)  has  finite  discontinuities 
at  certain  points. 

The  value  of  ,         2  cos  bDxf(x)  is  assumed  to  exist,  in  general, 

but  it  may  not  exist  at  certain  points.     It  need  not,  in  general, 
be  continuous  either  as  a  function  of  x  or  b. 

We  may,  therefore,  conclude  from  (35)  that:  Whenever  2  cos 
aDx  exists  and  is  a  continuous  function  of  x,  then  the  solution 
u(x)  is  a  continuous  function  of  x,  and  is  given  by  equation  (30). 
At  particular  values  of  x,  where  the  function  f{x)  or  its  ^^-derivatives 
of  any  order  have  a  finite  discontinuity,  2  cos  aDxf{x)  has  a 
finite  discontinuity,  and  this  is  also  the  discontinuity  of  the 
function  u{x}.  For  example,  suppose  that  in  the  neighborhood 
of  the  point  Xj,  the  function  f(x)  is  given  by  the  equation 

f{x)=(l){x)   when  x<Xi   and  f(x)  =(f)ix) +A+ Bix-Xj)* 


Snow]  Spectrum  Analysis  of  Radio  Signals  249 

when  x>Xi  where  (f>{x)  is  continuous  with  all  its  derivatives  at  x^. 
Then  f(x)  jumps  by  an  amount  A  as  we  pass  through  the  point  Xj 
and  its  fourth  derivativejumps  an  amount  \4.B.  Hence,  u(x)  jumps 
by  an  amount  2  (A  —a^B).  If  B  equals  zero,  the  jump  in  u{x)  is 
twice  that  in  /  {x) ,  and,  in  general,  the  discontinuity  in  u{x)  is  greater 
than  that  in  f{x),  although  in  this  case  if  B  =a^A  there  would  be  no 
discontinuity  in  u{x).  The  function  f(x)  would,  however,  be  con- 
tinuous if  A  =0,  together  with  its  first  three  derivatives,  and  yet 
u{x)  would  have  a  discontinuity  of  —2a^B.  The  ordinates  of 
the  curve  f{x) ,  its  slope  and  its  ciu-vature  would  all  be  continuous 
at  Xi  and  there  would  be  notliing  in  its  appearance  to  suggest  the 
discontinuity  in  the  spectrum. 

Example  illiistrating  this  method:  Suppose  for  example  that  the 
observed  curve  is  given  by 

/(x)=^[tan-^-tan-^] 

X2>Xi'>o  and  6>a 

/x  —  ^\ 

tan"^  I  -^T —  j  =imaginary  part  of  log  [6  -[-iix^  —  x)'] 


where 

Now 


Hence  if  a  <6 

2  cos  a  D^r  tan"^  /x„  —  x 


f^^\^(ei^D,  +e-iaD,)  tan-i  (\^) 

=imaginary  part  of  (e'^^x  +e~^'^^x)  log  [b  +i{x2  —  x)] 

=imaginary  part  of  log  [6 +*(^2~^~*'^)]  +  log  [b +i (X2  —  x +ia)] 

=imaginary  part  of  log  [b  +a  +i{x2—x)]  +log  [b  —  a  +i{x2  —  x)] 

=tan-i  ('^)  +  tan-i  f"^) 
\b+a/  \b  —  a/ 

Hence 

.  cos  a  D^  m  =^  [tan-  {^)-t.n-'  (^)] 
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This  gives  the  first  term  on  the  right  side  of  (30).     The  second 
term  on  the  right  side  of  (30)  is  the  integral. 

00 
Aa  C         dx'         r^     .,X2—x'     ^     _,  x.-acn 

—00 

It  is  next  necessary  to  evaluate  these  integrals.     To  do  this  note 
that  the  following  is  identically  true : 


^[--'e^')--e-^')]=^j 


ds 


b'+{s-xy 

Multiplying  this  identity  by ^     .  ,  _  .    and  integrating  from 

%'==  —  £»  to  a;'  =  +  00  gives  the  term  desired  as  the  double  integral 

00  12 

A  ah  C        dx'         C         ds 


\ah  C        dx'         C 

x^  J  a^-V{x'-xY]  ^^-{s-x') 

—00  X\ 


{h"  +{x'  -sy\\a-  ^{x'  -xy\ 
=  ^ii^±^  f ^£ ■  (by  34) 

Hence,  the  formal  solution  (30)  gives  the  spectral  distribution 
of  energy  u  {x)  by 

CO 

w  (x)  =  2  cos  a  D^f{x)  ~  -       2  ,  /  / \2 

'^        IT  J  a^  +{x'  —xy 

—  CO 

when  the  observed  function  is 

/  (^)  =  ^  ftan-^  (^X^)  ~  t^""'  (^T^)l  "^^^""^  ^ >  "^ 
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These  two  curves  are  of  the  same  type  and  are  somewhat  similar 
to  the  curve  showing  the  magnetic  field  along  the  axis  of  a  solenoid. 
The  observed  curve  is  broader  than  the  energy  curve.  The  latter 
becomes  infinitely  sharp  at  the  edges  x^  and  X2  when  b  approaches  a. 
In  this  case,  the  energy  distribution  is  zero  outside  the  range 
Xi<x<X2  and  has  the  constant  value  A  within  that  range.  The 
observed  curve,  however,  approaches  a  limiting  curve  which  is 
continuous.  The  parameter  b  —  a  which  occtus  in  the  energy 
curve  is  vaguely  analogous  to  what  is  called  the  decrement  of  the 
wave.  This  must  be  increased  by  the  damping  coefficient,  a,  of 
the  receiving  circuit  in  order  to  obtain  tlie  observed  curve.  This 
enlargement  of  what  is  usually  called  decrement  is  dealt  with  in 
more  general  fashion  later.  In  many  simple  types  of  waves  it 
represents  a  quantity  in  the  wave,  which  is  augmented  by  a 
when  received  by  the  circuit. 

This  third  form  of  solution  (30)  like  the  two  previous  forms  has 
the  advantage  that  it  may  be  applied  to  cases  where  the  function 
f{x)  is  represented  by  different  analytical  expressions,  in  different 
parts  of  its  range,  which  may  or  may  not  fit  on  to  each  other 
smoothly.  It  involves  graphical  differentiation  and  integration. 
To  apply  the  second  form  of  solution  requires  harmonic  analysis 
which  may  be  affected  mechanically.  However,  it  is  sometimes 
the  case  that  less  labor  is  involved  by  representing  the  observed 
function  by  an  empirical  equation  so  as  to  gain  the  advantages 
which  analytical  tools  often  possess  over  mechanical  ones.  In 
order  to  find  what  types  of  elementary  solutions  might  be  most 
appropriate  for  empirical  representation  of  the  observed  function, 
a  few  simple  types  of  elementary  solutions  will  be  considered. 

V.  FOURTH   FORM    OF   SOLUTION— EMPIRICAL   FORMULA 
FOR   THE   OBSERVED   FUNCTION. 

A  very  simple  expression  for  the  energy  distribution  results 
when  the  observed  function  /(x)  is  of  the  simple  type 


n  a\  +iXn-Xy 


where  an>a  and  x^  is  any  constant.  The  solution  u(x)  is  given 
by  simply  replacing  an  by  the  positive  quantity  a^  —  a  in  this 
equation,  a  form  which  is  generalized  later.  This  result  follows 
as  a  special  case  of  the  solution  (9) .     If  we  place 

'  a\  +{Xn-x'y 
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in  (9)  it  gives 

CO  00 

u(x)  = I  -^ — -7-7 r;  I  e^"  COS  «  (x'  —x)da 

TT    J  a\+{x'-x^Yj 

—00  0 

where  the  integration  must  be  performed  first  with  respect  to  x' 
in  order  that  this  may  have  a  meaning.  To  do  this  let  x'  =Xt,  + 1. 
Then 


«w-^J...„J?^ 


-x+t)dt 


+  P 

00 
"cos  ex  t 


= \  e^"  d  a  cos  a  (Xn-x)  I  - 

TT      J  Jan- 

0  —00 


dt 


sin  at  ,, 
dt 


A  an  C  ^„  ,       .        /  \  fsi 

0  —00 

00  00 

T,   .  r  sin  «  i  ,^     „       J  fcos  a  ^  j^     tt 

But  -TT^dt  =0  and       „    ,  ,,  dt  =  —  e'^" 

—CO  —00 

(Laplace's  integral).     Hence 

OD 

u{x)  =AS^"-  d  a  e'^"  cos  a  (xn  —  x) 


=  Ay*e-(^-'-^^"  cos oi  (xn-x)  da 
0 

A  (an -a) 

....  if  an  >  a 


{an-ay  +ixn-xy 
=  meaningless  if  an  <  a 


It  is  easy  to  see  that  if  a  pair  of  fmictions  f(x)  and  u(x)  are 
known,  such  as 

f(x)  = "^^  ^° 

^^""^     a'+(xn-xr 
and  (36) 

AnCOn-o) 

(an-a)-+(xn-^)^ 

which  satisfy  the  original  equation  (5),  then  other  pairs  of 
functions  which  will  satisfy  (5)  may  be  obtained  by  differentiating 
this  pair  with  respect  to  an,  Xn,  or  x ;  or  by  integrating  with  respect 
to  these  quantities.  Also,  by  making  any  linear  combinations  of 
these  simple  solutions,  other  simple  solutions  may  be  found. 


Stww] 
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The  simple  type  of  solution  here  given  is,  perhaps,  the  one  of 
most  practical  value  for  empirical  representation  of  the  observed 
function. 

It  might  be  noticed  that  this  type  of  solution  is  not  strictly 
compatible  with  the  fact  previously  noted  that  f(x)  and  u{x) 
are  both  odd  functions  of  x,  whereas  the  pair  just  exhibited  are 
neither  even  nor  odd  functions.  By  changing  the  sign  of  Xu 
and  subtracting,  however,  we  get  an  odd  function 


=  4Xn  An 


anX 


=  4JCn  A 


[X^  +Xn  +a\  -  2XnX]  [X^  -hX^ri  +a\  +2XnX] 

anX 

°  /'-v2 


(x^  +Xa  +a\)  ^  -  4xhj,' 


which  is  an  odd  function  of  x,  and,  therefore,  corresponds  to  a 
physical  solution.  However,  as  we  have  seen,  the  extra  term 
here  added,  namely, 

I 

a^  +(xn+xy 

is  utterly  negligible  compared  to  the  term 

I 

a*  +  {Xn  -  x)  ^ 

when  we  are  concerned  only  with  a  narrow  range  of  high  positive 
values  of  the  frequency,  and,  therefore,  the  simpler  type  of  solution 
is  the  more  practicable. 

By  differentiation  and  integration  of  the  fundamental  types 
(36)  many  other  simple  types  are  obtainable,  such  as 


I/W  = 


A  an 


(x-Xny+a\ 
II/W=A[tan-.(S^)  +  tan-.(^)] 

A{x-Xn) 


III  fix) 


{x-XnY  +a'^n 

IV  fix)  =A  logT^^^Il^^i^^+^l 
'^^  \ix-x,)'+a\j 


(37) 
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where 

an>a 
or  more  generally 

t\x-x')A{x')-^a^B{x') 

where  A  and  B  are  any  functions  of  x'.  Other  forms  involving 
any  power  of  the  denominator  may  also  be  obtained.  The  corre- 
sponding solution  u{x)  for  all  of  the  above  types  (or  of  sums  and 
Hnear  combinations  of  them)  is  found  by  substituting  a^  —  a  for 
fln  in  these  expressions  for  f{x). 

This  rule  is  a  generalization  of  the  principle  which  forms  the 
basis  of  all  decrement  measurements.  It  corresponds  to  the  fact 
that  peaks  in  the  curve  f{x)  correspond  to  peaks  in  the  u{x)  curve, 
which  are  narrower  by  2a,  the  breadth  of  the  resonance  peak  of 
the  receiving  circuit. 

The  physical  meaning  of  these  types  may  be  illustrated  by  con- 
sidering the  radiation  produced  by  another  simple  circuit  at  a 
distance,  oscillating  freely  and  sending  out  a  radiation  which  is 
given  by  (S  {t)  =Ce~"^^  cos  2'n:x't.  The  development  of  this 
function  of  the  time  in  a  Fourier  integral  for  all  values  of  the  time 
from  —  00  to  +  00  shows  that  the  amplitude  of  each  harmonic  com- 
ponent of  frequency  x  is  given  by 

C  x'  +x"  +b' 

^^^     2-K'[{x-x'y+b'][{x+x'y-Vh''] 

very  approximately 


where 


[(x-x'y+b'] 

c 


B  = 


47r*b 
Since  u(x)^x  E^  {x)  this  gives 

Bbx  Bbx 


ix-x'Y+b''  '  {x+x'y+b'' 

Bb(x-x')  +x'Bb     Bb{x+x')-x'Bb 

{x-x'Y+b'      ^       {x+x'Y+b^ 

A  {x-x')  -Vx'Bb     A  {x+x')-x'Bb 


{x-x'y+¥      '       {x+x'y+b^ 
where  A  =Bb 
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According  to  the  rule  above  given,  the  observed  function,  when 
the  radiation  for  the  simple  circuit  is  received  by  one  whose  damp- 
ing constant  is  a,  would  then  be  given  by 

fix)  -8RLx  ^Hx)     -4(x-xO  +x'B(b+a)     A(x+x') -x' B{b +a) 

exactly 

A(x-x')+x^B(b+a)  _^    (b_+a^+4^^ 
{x-x'y+b+ay       -^^  ix-x')'+{b+ay 

approximately 

=x'B  -, ,.,  ,^,    ,    ,„   if  the  range        ,^^  is  small. 

{x  —  x)^+{b+a)^  x 

This  is  a  combination  of  the  types  (37),  I  and  III.  Therefore 
the  oscillations  of  a  distant  simple  circuit  which  produces  an 
effective  current  in  the  receiving  circuit  of  the  form 

^  ^  ^     SRLx    {x-x')^+a\ 
or  approximately 

B  ttn 


8RL^^  {x)  =i 


'[ix-xy+a\] 

when  the  range  -^ ~  is  small,  must  have  a  spectral  energy  dis- 
tribution 

8  RLE'  (x)  =7 ^xfrr""^    ^. 

^  ^      ix-x'y+{a^-ay 

The  width  of  the  peak  of  this  curve  is  2(an  — a);  that  is,  it  is 
narrower  than  the  peak  of  the  observed  current  by  2a,  the  width 
of  the  resonance  peak  of  the  receiving  circuit.  The  receiving  cir- 
cuit adds  its  own  decrement  to  that  of  the  radiation. 

The  symmetrical  simple  solution  of  the  type  (37),  I,  is  therefore 
seen  to  correspond  to  one  of  the  most  important  physical  types 
of  radiation  to  be  met  with,  namely,  damped  oscillations.  It  is 
readily  seen  to  contain  the  case  of  undamped  oscillations  as  a 
limiting  case,  for  if  On  is  but  very  slightly  greater  than  a,  the 
solution 

a^  —  a 
{x-x^y  +{ar,-ay 
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is  very  small  for  all  values  of  x,  except  x  =Xa,  which  corresponds 
to  a  single  spectrum  line  of  frequency  Xn',  that  is,  an  undamped 
wave.  The  unsymmetrical  form  (37),  III,  which  changes  sign  at 
X  =Xn  is  also  seen  to  possess  a  similar  property.  It  approaches 
zero  in  general,  but  becomes  +  or—  00  when  x  approaches  Xa  from 
opposite  sides,  when  an  approaches  a  as  a  limit.  Although  simple 
solutions  which  change  sign  do  not  correspond  to  any  physical 
case  when  considered  alone,  yet  a  number  of  them  together  may 
readily  be  compounded  into  a  function  which  is  everywhere  posi- 
tive inside  a  given  range,  just  as  in  the  case  of  sines  and  cosines. 
This  type  of  solution  is  therefore  useful  in  representing  the  ob- 
served function.  The  form  (37),  IV,  also  changes  sign  half  way 
between  x^  and  X2.  The  form  (37),  II,  is  very  useful  in  some 
cases.  It  is  easily  seen  that  if  a^  approaches  a  as  a  limit,  the 
solution  u{x)  is  zero  everywhere  outside  the  range  between  x^  and 
X2,  and  is  constant  inside  this  range. 

VI.  GENERAL  RELATION  BY  MEANS   OF   COMPLEX 
INTEGRATION. 

All  of  these  special  cases  are  included  in  a  very  general  relation 
which  may  be  derived  from  (30) ,  and  which  gives  an  elegant  solu- 
tion in  case  the  function  f(x)  has  been  represented  by  an  analytical 
expression.  If  the  poles  of  f(z)  all  lie  at  a  distance  from  the  real 
axis  (in  the  complex  2;-plane),  which  is  greater  than  a,  then  2  cos 
aD^  fix)  or  {e^^^^  +e~^^^)  fix)  represents  fix+ia)  +fix  —  ia)  since 
the  expression  e^^^^^  fix)  is  just  the  symbolic  form  of  Taylor's 
theorem,  and  in  this  case  the  expansion  is  valid  about  the  point  x, 
which  is  any  point  on  the  real  axis.     Equation  (32)  is  therefore 


cx> 

u  {x)  =f  ix  +ia)  +f  ix  —  ia) I  — 


/  jx')  dx' 

+ix'-xy 


If  f^iz)  represents  that  part  of  fiz)  whose  poles  are  in  the  upper 
half  plane  and  /jC^)  that  part  whose  poles  are  below  the  real  axis, 
and  if  fiz)  is  finite  or  zero  at  infinity,  then  by  deforming  the  path 
of  integration  indefinitely  upward  in  the  one  case  and  downward 
in  the  other  case,  and  evaluating  the  integrals  around  the  poles 
of  the  nucleus  by  Cauchy's  theorem  of  residues,  one  easily  arrives 
at  a  fourth  form  of  solution 

«  (*)  =/i  (^  +ia)  +/2  (^  -  i(^)  (38) 
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This  means  that  if  a  new  function  is  formed  from  f(x)  by  placing 
all  its  poles  nearer  to  the  real  axis  by  an  amount  a,  the  new  func- 
tion is  u{x).     B}^  similarly  treating  the  original  equation  (5)  one 

finds 

f{x)  =u-^{x  —  ia) —u^ix -\-ia) 

which  means  that  if  a  new  fimction  is  formed  from  u{x)  by 
moving  all  its  poles  a  distance  a  farther  from  the  real  axis,  this 
new  function  is  fix) . 

VII.  APPLICATION    OF   FOURTH   METHOD   TO 
EXPERIMENTAL    CURVES. 

In  the  analysis  of  a  given  radiation  one  is  usually  only  concerned 
with  relative  values  of  the  induced  current.  When  the  frequency 
range  studied  is  not  large  compared  to  the  mean  frequency,  one 
may  obtain  correct  relative  values  by  identifying  fix)  with  ^^{x) 
and  u{x)  with  E-ix),  for  multiplying  ^^{x)  by  x  distorts  the  curve 
in  a  manner  which  is  practically  all  compensated  for  when  we 
divide  u{x)  by  x  to  obtain  E^"ix). 

When  the  capacity  of  the  circuit  is  varied,  a  wave  meter  may  be 

used  to  find  the  corresponding  resonance  frequency  X,  which  is 

practically  the  same  as  the  characteristic  frequency  x  =  -y/X^  —  a^. 

This  method  then  gives  the  true  frequency  which  depends  not 

only  upon  the  capacity  of  the  condenser,  but  upon  the  inherent 

distributed  capacity  of  the  coil. 

R  X 

The  constant  a  = =-  (  = —  5  where  5  is  the  decrement  of 

^  IT  L  27r 

the  circuit)  must  be  determined  at  a  value  of  the  frequency  near 
which  it  is  to  be  used.  Since  the  resistance  and  inductance  of  the 
circuit  will  change  with  frequency,  due  to  eddy  currents,  the  value 
of  a  will  not  be  strictly  constant,  but  may  be  considered  so  for 
relatively  small  ranges  of  frequency.  If  an  unmodulated  (simple 
harmonic)  wave  of  frequency  x'  be  studied  with  the  circuit,  the 
current  curve  should  be  of  the  form 


so  that 


If  {x  —  x')"^  be  plotted  against  ^^  .  .>  and  the  result  is  a  straight 


3^^ 

(x)  = 

{X- 

xy+a' 

-xy 

+a^ 
I 

A 

^Hx) 
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line,  the  intercept  of  this  line  gives  the  value  of  a^  to  be  used  in  the 
neighborhood  of  the  frequency  x'. 

Figure  i  illustrates  the  use  of  an  empirical  formula  of  the  simple 
type: 

where 

A(x)=A.[tan-^(^-^^)  +  tan-^(^)] 


fM- 


^2     J      \(^~^^^ 


A, 


J{X)=X^{71)  IN  ARBITRARY  UNITS  (OBSERVED  FUNCTION) 


moo 


1900 


2000  — -X 

X'FREQUENCV  IN  KILOCYCLES  PER  SECOND 


Fig.  I . — Energy  distribvtion  in  a  condenser  modulated  wave. 

The  observed  values  are  represented  by  circles  where,  as  in  the 
previous  notation,  f.{x)  is  x  3^  W,  and  3^  is  the  effective  ciurent 
produced  in  the  circuit,  as  indicated  by  a  hot-wire  ammeter. 
The  imits  of  current  are  entirely  arbitrary,  since  only  relative 
values  are  desired.  The  radiation  in  this  case  was  a  "condenser 
modulated,"  one  caused  by  the  rotation  of  the  plates  of  an  air 
condenser  in  a  sending  circuit  in  which  there  is  a  sustained  simple 
harmonic  wave.  The  experiments  were  performed  by  R.  T.  Cox 
and  Miss  F.  Kenyon,  of  the  Bureau  of  Standards  radio  laboratory. 
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The  numerical  value  of  the  constants  which  were  used  to  fit  the 
observed  curve  were: 
Ai  =0.672(10)',  x\  =1857.5  kilocycles  per  sec.  a^  =12.5  10'  sec.~^ 

Xi  =1982.5  kilocycles  per  sec. 
A3  =2.56(10)^      X2  =1882.0  kilocycles  per  sec.  Oj  =i6  10^  sec."^ 
A3  =  —  1.10(10)",  X3  =1948.0  kilocycles  per  sec.  a^  =22  10*  sec.~^ 

The  circuit  constant  a  =10.7  10^  sec."^  was  determined,  as 
explained  above,  from  the  curve  obtained  on  the  unmodulated 
wave.     The  inductance  L  was  56.8  microhenrys. 

The  empirical  formula  in  this  case  fits  the  observed  curve  with 
no  more  error  than  those  of  experiment. 

The  energy  distribution  was  computed  from  u(x)  where  u{x) 
was  obtained  by  decreasing  each  a^  by  a  in  the  empirical  formula 
for  f(x).  It  is  represented  by  the  dotted  ciu*ve,  and  shows  that 
the  spectrum  of  the  radiation  consisted  of  a  wide  band  of  nearly 
uniform  intensity. 

The  representation  of  a  function  which  has  well-defined  peaks 
by  an  empirical  equation  of  the  type 

2^n  I 


is  generally  fairly  easy  provided  that  each  an  is  greater  than  a. 
If  the  peaks  are  no  wider  than  2a,  this  would  indicate  the  presence 
of  a  line  spectrum.  This  was  the  conclusion  drawn  from  an 
attempt  to  represent  the  curves  of  Figure  2  by  terms  of  the  above 
type.  The  two  curves  are  drawn  through  the  observed  values  for 
two  different  values  of  the  resistance  of  the  receiving  circuit,  the 
radiation  being  the  same  in  both  cases.  It  was  produced  by  inter- 
rupting a  simple  harmonic  wave  with  a  chopper  whose  frequency 
was  612  cycles  per  second.  The  frequency  of  the  wave  which  was 
modulated  was  53  kilocycles  per  second.  One  would  expect  that, 
if  the  amplitude  of  a  pure  harmonic  vibration  were  varied  harmoni- 
cally with  a  frequency  of  612  cycles  per  second,  the  frequency  of 
the  original  note  would  be  absent  from  the  spectrum,  and  there 
would  appear  frequencies  of  this  value  plus  and  minus  multiples 
of  the  612,  so  that  the  spectrum  would  consist  of  lines  separated 
by  1,224  cycles.  The  separation  of  the  lines  is  about  1,200 
cycles,  but  the  amplitudes  are  not  symmetrical,  showing  that  the 
action  of  the  chopper  was  not  so  simple  as  this. 
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The  crosses  in  Figure  2  represent  the  computed  currents  on  the 
assumption  that  the  spectrum  consisted  of  three  monochromatic 


9  (X)  =  E.FrECTIVE  CURRENT 


53  54  55  56  57 

X=  FREQUENCY  IN  KILOCYCLES  PER  SECOND 

Fig.  2.— Current  in  two  different  circuits  produced  by  the  same  u-ave. 

lines  at  52.5,  53.7,  and  54.9  kilocycles  per  second,  respectively 
The  equation  used  was 

I 


K=3     A 


„_,    a   ^   ^  /x-XnV 


Where 


a  =0.206  10  ^  in  lower  curve 
a  =0.65   10  ^  in  upper  curve 


The  agreement  of  computed  values  with  observ^ed  values  in  the 
case  of  the  circuit  which  had  the  higher  resolving  power  justifies 
the  inference  as  to  a  line  spectrum.  The  deviations  of  computed 
values  from  observed  values  in  the  case  of  the  upper  curve  might 
be  accounted  for  by  a  slight  underestimate  of  the  resistance  of 
this  circuit.  The  resistances  in  both  cases  were  determined  from 
decrement  measurements.  All  three  curves  were  obtained  by  the 
same  observ^ers  in  a  preliminary  study,  the  object  of  which  was  the 
development  of  a  method  by  which  the  interference-producing 
quality  of  a  wave  might  be  characterized. 
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VIII.  CONCLUSION. 

Four  different  methods  have  been  shown  for  finding  the  energy 
spectrum  of  a  radiation  in  terms  of  the  mean  current  it  produces 
in  a  receiving  circuit,  this  current  being  observed  for  all  real 
values  of  the  characteristic  frequency  of  the  circuit.  The  problem 
consists  essentially  in  finding  practicable  methods  of  solving  the 
integral  equation  for  u{x) 

CD 

m(x')  dx' 


m-ll 


—  (DO 

when  /(%)  is  given  as  an  odd  function  of  x. 

It  is  obvious  that  other  applications  of  resonance  theory  may 
involve  solving  this  equation,  which  may  be^  quite  different  from 
the  purely  electrical  problem  with  which  it  is  here  associated. 

The  first  method  gives  the  solution  of  this  equation  in  terms  of 
a  double  integral  analogous  to  Fourier's  (equation  10).  This 
form  is  only  of  value  for  computing  when  /(x)  is  given  analytically, 
but  is  useful  as  a  basis  for  deriving  other  forms. 

The  second  method  is  based  upon  the  fact  that  when  the  nucleus 
of  an  integral  equation  is  an  even  function  of  the  sum  or  difference 
of  the  two  variables  its  corresponding  normal  functions  are 
trigonometric  functions.  Harmonic  analysis  of  the  nucleus 
(mechanically  or  graphically)  amounts  to  a  determination  of  the 
characteristic  constants.  A  similar  analysis  of  the  observed 
function  has  been  shown  to  be  sufficient  to  lead  to  a  series  solution 
(equation  26). 

The  third  form  gives  the  solution  in  terms  of  differentiations 
and  integrations  of  the  given  function  /(x).  Like  the  first  two 
forms  of  solution,  it  is  applicable  to  cases  where  the  functions 
fix)  or  u{x)  possess  finite  discontinuities,  and  has  the  additional 
advantage  that  the  operations  may  be  performed  graphically 
(equation  30). 

The  fourth  method  consists  in  representing  the  observed 
function  by  empirical  equations  of  certain  simple  types  (37),  the 
solution  corresponding  being  easily  obtained  by  subtracting  the 
damping  coefficient  of  the  circuit  from  certain  analogous  para- 
meters in  the  simple  types  of  solutions  used.  This  m.ethod  is 
applicable  to  the  limiting  case  of  line  spectra.  It  is  applied  to 
three  different  experimental  curves. 

Washington,  July  3,  1923. 
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